Heavy Fermion Quantum Criticality 
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During the last few years, investigations of Rare-Earth materials háve made clear that not only 
the heavy fermion phase in these systems provides interesting physies, but the quantum criticality 
where such a phase dies exhibits novel phase transition physies not fuUy understood. Moreover, 
attempts to study the critical point numerically face the infamous fermion sign problém, which 
limits their accuracy. Effective action techniques and Callan-Symanzik equations háve been very 
popular in high energy physies, where they enjoy a good record of success. Yet, they háve been 
httle exploited for fermionic systems in condensed matter physies. In this work, we apply the RG 
effective action and Callan-Symanzik techiques to the heavy fermion problém. We write for the 
first time the effective action describing the low energy physies of the systém. The /-fermions 
are replaced by a dynamical scalar field whose nonzero expected value corresponds to the heavy 
fermion phase. This removes the fermion sign problém, making the effective action amenable to 
numerical studies as the effective theory is bosonic. Renormalization group studies of the effective 
action can be performed to extract approximations to nonperturbative effects at the transition. By 
performing one-loop renormalizations, resummed via Callan-Symanzik methods, we describe the 
heavy fermion criticality and predict the heavy fermion critical dynamical susceptibility and critical 
specific heat. The specific heat coefiicient exponent we obtain (0.39) is in excellent agreement with 
the experimental result at low temperatures (0.4). 

PACS numbers: 71.27.-|-a, 73.43.Nq, 71.10.Hf 
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For a couple of decades, heavy fermion materials háve 
attracted the focus of a large part of the experimental and 
theoretical condensed matter commumty[l|, [4]. There 
are many reasons for such a spothght on these materi- 
als. They exhibit exotic superconductivity, interesting 
magnetism, but most importantly heavy quasiparticles 
with an enlarged Fermi surface. This heavy quasipar- 
ticle phase perishes into a quantum critical point with 
interesting, puzzling and not yet understood nature[3|,|j]. 

It is quite striking to see the fermionic quasiparticle 
with masses from tenths to about thousandths of an el- 
ementary electron mass. This has been understood as 
arising from Kondo-like physies of the almost localized 
/— electrons when they hybridize with the lower atomic 
angular momentům bands of the materiály, |5|. This hy- 
bridization gives rise to an enlarged Fermi surface as the 
/—electrons now contribute to the Fermi volume, and to 
the large quasiparticle mass and large specific heat coeffi- 
cient, through enhanced coUective Kondo-like low energy 
scattering. Of course, the /—electrons háve a strong ten- 
dency to localize due to their large U which fights the 



hybridization V . 

Some of these materials can be tuned (by apply- 
ing pressure, etc.) so that a critical value {U /V)c is 
reached such that, for values larger than the critical 
value, the /—electrons localize and there is no heavy 
fermion phase. Instead, there is a small Fermi surface 
metal that usually exhibits magnetic order mediated via 
RKKY interactions'6'|. For subcritical values of U /V , the 
systém is in the heavy fermion phase with large Fermi 
surface and no magnetism. 

At the critical value, a continous quantum phase tran- 
sition occurs as corroborated via scaling experiments. 
This transition is not understood. The lack of under- 
standing is a barrier to the fuU characterization of the 
physical properties, phase diagram and experimental re- 
sponse features of these materials. In this work, we turn 
our attention to understanding this heavy fermion quan- 
tum criticality ji, 01 . 

We start from the partition function for /—electron 
hybridizing with metallic c— electrons, i.e. the periodic 
Anderson Model 
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where the subscripts M, L (and all capital letters sub- 
scripts for that matter) indicate the angular momentům 
degeneracy of the ground statě (j = 5/2), U is the Hub- 
bard repulsion and V represents the strength of mixing 
of /-electrons with the conduction band. 

If we write n^ = ^^ ÍJmÍím, the Hubbard interaction 
term then takés the form (in Euchdean time and reál 
space) ^ ^j / dr [n? — Ui] . We can decouple this 
interaction term by using the Hubbard-Stratonovich 
identity g-Z-^^S^f"? = Jp<^e- Z"*"^^ ^[^^+''^'"'] 
which in Euchdean time and momentům space 

Using 
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In order to háve / levels occupied we choose Eq < fi. 

This last action is what quantum field theorists would 
call the bare action. It embodies the essential physics 
of the heavy fermion phase and the heavy fermion crit- 
icahty. Unfortunately, it is very hard to solve the bare 
action exactly or numerically to high accuracy in order 
to extract the desired information from it. The famous 
fermion sign problém thwarts numerics, and exact Solu- 
tions are normally impossible in many body problems 
just as this one. Fortunately, there is a way forward that 
can help extract some, and perhaps a lot of the physics. 



The renormalization group|8|, and in particular effective 
action and Callan-Symanzik techniques[9l. IlOl. Ill| . popu- 
lar in particle physics, provide room for progress. 

One of the lessons of the renormalization group is that 
as we concentrate on longer wavelength, lower energy de- 
grees of freedom, the short distance and high energy fluc- 
tuations do two things. These fluctuations renormalize 
the strength of the terms in the originál action and they 
generate new terms in the action which in turn change as 
they get renormalized. The end result is that some terms 
in the action bccome larger while others become smaller, 
thus not contributing to the to the universal low energy 
physics of the systém. We will thus analýze the action for 
heavy fermion materials above, and obtain the effective 
action with terms relevant to the low energy physics of 
the heavy fermion phase and to the critical point where 
such a transition perishes. 

To obtain the effective action, rather than obtain the 
renormalizations all at once, it proves advantageous to 
integrate out the /—electrons and get those terms they 
contribute to the effective action which are relevant for 
the low energy universal physics of the critical point and 
heavy fermion phase. This is different from the tradi- 
tional approach where the /—electrons are integrated out 
and all terms are kept, usually in the form of a functional 
determinant, irrespective of whether the terms are rele- 
vant or not for the low energy physics. The other im- 
portant novelty in our analysis is the application of RG 
effective action and Callan-Symanzik mcthods to access 
the heavy fermion critical point from the effective action. 

After the calculations are performed, the /—electrons 
disappear from the theory. We háve instead the 
Hubbard-Stratonovih field (/?, which has acquired dy- 
namics through the /— electron fluctuations. Such 
fluctuations also generate self-interaction terms for the 
íy3-field, and interaction terms between the ip^s and the 
metallic c— electrons. 

The effective action for the heavy fermion materials 
comes out to be 
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We see that the universal physics of the heavy fermion 
systém is captured by an action of dynamical scalar fiels 
interacting with the metallic c— electrons. We call this 
action the heavy-fermion dynamic ip'^ action. The heavy 
fermion phase corresponds to 32 being negative and ip 
acquiring a nonzero expected value as ip is proportional 
to the density of /—electrons that hybridizes with the 
metahic ones. The heavy fermion critical point occurs at 
f/2 = O, when {(p) first becomes 0. 

The heavy-fermion dynamic </?* action is a new and im- 
portant result. It opens the dooř to accurate numerics for 



the fermion action, as the interacting fermions that drive 
the transtion háve been replaced by a scalar field. This 
should eliminate the fermion sign problém that plagues 
numerics, for all the action is happening in the scalar 
fields and not the left-over metallic fermions. 

One can apply standard order parameter RG to this ac- 
tion. As an example, below we do a one-loop momentům 
shell renormalization with the help of Callan-Symanzik 
equations á la Weinberg ll|, [l^, ll3(to resum and thus 
catch some of the nonperturbative physics of the transi- 
tion. This can of course be improved by going to higher 
orders, and there is also plenty of room to perform e- 
expansion studies instead of momentům shell. 



Since we can use the bare Fermi velocity of the metal- 
lic c— electrons as a standard of speed in the materiál, we 
use it as such to express our frequcncies in units of mo- 
mentům and work in "God-given" heavy fermion units: 
vp = 1 and h = 1, kp — m — A. After renormalization, 
the heavy-fermion dynamic ip'^ action becomes 
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We now move to consider the specific momentům shell 
renormalizations to determine the renormalization fac- 
tors. The inverse p propagátor, uj — Ag2 goes into 
-IJ,g2 Z^ + -^^w -I- Zkkr = -A52 + t^ + S^, where 
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Next we tackle the renormalization of the c— elec- 
tron inverse propagátor w — (e^: — /i) which becomes 

Zl [lo - (e^ - ^)] =uj - (e^ - A^) + Sc(w, k), with 
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Notice that the c— electron propagátor only renormalizes 
trivially to this order at the critical point. The reason 
is that if (72 is nonzero, we need to include the nonzero 
expected value of ip and expand around it in the heavy 
fermion phase in order to get the physics right. We are 
not interested in the heavy phase but in the critical point, 
where such a phase disappears. At the critical point, 
g2 — O and things work out as expected. 

We now turn our attention to vertex renormalizations, 
starting with the p — c vertex, A^'^ g^. It renormalizes 
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And finally to the renormalization of the Lp interaction, 



fj,-^ gf Z't = A-254 + r^, with 



Such a propagátor is proportional to (1^9) , which satisfies 
a Callan-Symanzik equation ll|, [l^] 
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The heavy fermion criticahty occurs when 52 = 0. In 
the 52 = O critical manifold 
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The iy9 anomalous dimension at cricitcahty is 
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We can make experimental predictions of our heavy 
fermion ip^ action. These predictions are part of our new 
results and serve as an experimental test of the validity of 
our theory. Among those are the susceptibility and spe- 
cific heat of the systém. In the heavy fermion phase, the 
/-electrons contribute to the susceptibility and specific 
heat since they are mixed with the metallic band elec- 
trons, forming in conjunction one "Fermi liquid" . The 
susceptibility of conduction electrons is negligible since 
they are nonmagnetic. The susceptibility is proportional 
to the number of /-electrons being "pulled" by an applied 
field, and hence to their density of states, which is given 
by the imaginary part of the /— electron propagátory]. 



The specific heat coefficient is also proportional to the 
density of states, which is inversely proportional to the 
Fermi energy[2]. The /-electrons, being quasi-localized, 
form a quite thin band and hence háve a small Ep ■ Thus 
their density of states is so big in comparison with that 
of the conduction electrons, that the density of states of 
the "mixed Fermi liquid" can be approximated to be that 
of the /-electrons. We obtain 
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We obtain a specific heat coefficient exponent 7i^/2 = 
0.39. This is in excellent agreement with the exponent 
0.4 found for YbRh2SÍ2 at low temper atures[3|. 

Via renormalization group studies and effective ac- 
tion techniques common to field theories of particle 
physics, we háve obtained the effective field theory for 
heavy fermion quantum criticahty. This marks impor- 
tant progress as the effective field theory is bosonic, viti- 
ating the fermion sign problém and thus being amenable 
to numerical studies and high order e expansion studies. 

The critical field theory can be studied using the renor- 
malization group. We did so via one-loop renormalization 
studies, improved by means of Callan-Symanzik resuma- 
tions to access some of the nonpertrubative effects. We 
thus make predictions for the exponents that character- 
ize the critical divergence of the specific heat coefficient 
and the critical charge susceptibility. Our specific heat 
coefficient exponent of 0.39 is in excellent agreement with 
the 0.4 found in experiments at low temperatures[3|. 
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